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Abstract 

A comprehensive analysis of general relativistic spacetimes which admit a shear-free, irrotational and 
geodesic timelike congruence is presented. The equations governing the models for a general energy- 
momentum tensor are written down. Coordinates in which the metric of such spacetimes takes on a 
simplified form are established. The general subcases of 'zero anisotropic stress', 'zero heat flux vector' 
and 'two component fluids' are investigated. In particular, perfect fluid Friedmann-Robertson- Walker mod- 
els and spatially homogeneous models are discussed. Models with a variety of physically relevant energy- 
momentum tensors are considered. Anisotropic fluid models and viscous fluid models with heat conduction 
are examined. Also, models with a perfect fluid plus a magnetic field or with pure radiation, and models 
with two non-collinear perfect fluids (satisfying a variety of physical conditions) are investigated. In partic- 
ular, models with a (single) perfect fluid which is tilting with respect to the shear-free, vorticity-free and 
acceleration-free timelike congruence are discussed. 
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1. Introduction 



In this paper, we shall consider general relativistic spacetimes which admit a shear-free, irrotational and 
geodesic (SIG) timelike congruence. It is well known pi ^ that if the stress-energy tensor is a perfect 
fluid whose flow lines form a SIG timelike congruence and whose density and pressure, as measured by a 
comoving observer, satisfy a barotropic equation of state, p = (with /i + p =/= 0), then the spacetime 
must be a Friedmann-Robertson- Walker (FRW) model. FRW spacetimes admit a six-parameter isometry 
group, Ge, whose orbits form spacelike hypersurfaces of constant curvature; coordinates can be chosen in 
such spacetimes so that the metric may be written as: 



ds 2 = -dt 2 



H 2 (t) 



dr 2 



(d9 2 



''9d4> 2 ) 



(1 



r2)2 



(1.1) 



where the constant n has been normalised to take the values ±1 and 0. It is also known [jlj, |3|] that if an 
equation of state p ~ p([i) is not assumed a priori then the resulting spacetime is still FRW. However, in 
the case of more general stress-energy tensors, there do exist solutions of the Einstein field equations that 
are not necessarily FRW models but where the spacetime admits a SIG timelike congruence. 

It is our aim here to explore the richness of spacetimes admitting a SIG timelike congruence. 
For illustration, we begin by presenting an example of a non-FRW spacetime that contains a SIG timelike 
congruence. Consider the following Bianchi VIo spacetime with "equal scale factors" (spatially homogeneous 
spacetimes will be studied in detail later): 



ds 2 = -dt 2 



X 2 (t)(dx 2 



dy 2 



dz 2 ) 



(1.2) 



Taking the source to be a fluid with four- velocity u a =5%, we observe that the fluid flow lines form a SIG 
timelike congruence. However, the source cannot be a perfect fluid since the Bianchi VIo models do not 
contain FRW models as special cases. Indeed, this model can be interpreted as an anisotropic fluid with 



/' 

P\\ 

P± 



3X 2 
~X 2 



1 

X 2 



2X 


X 2 


1 


~Y ' 


h X 2 


+ ^2 


2X 


X 2 


1 


~Y ' 


h X 2 


~ X 2 



(1.3) 
(1.4) 

(1.5) 



where /i is the energy density, and pii and p± are, respectively, the pressures parallel to and perpendicular to 
n a = XSq- Hence, there exists an anisotropic fluid solution whose flow lines are SIG. This example refutes 
the belief that spatially homogeneous spacetimes admitting a SIG timelike congruence must necessarily be 
FRW models. 

We are primarily interested in cosmological models. In such models (and, in fact, in all models) other 
physical conditions must be satisfied. The precise conditions depend on the particular energy-momentum 
tensor considered; for example, for viscous fluid models with heat conduction a suitable set of thermody- 
namical laws must also be satisfied. All models should satisfy the (various) energy conditions |||. We note 
that there exist models in the above example in which the weak, dominant and strong energy conditions 
are all satisfied. Consequently, physical constraints (or, at least, the energy conditions alone) will not force 
the models to be FRW. 

The remainder of the paper is organised as follows. In the next section the governing equations 
are displayed. For a spacetime in which there exists a shear-free, irrotational, geodesic timelike congruence, 
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coordinates can be chosen so that the metric takes on a simplified form [cf. equations ( [2.19 ) and (2.20)]; the 
governing equations are then given in these coordinates. In the following two sections general spacetimes 
are investigated with zero anisotropic stress and zero heat flux. At the end of section three, we show that in 
the special case of a perfect fluid whose flow forms a SIG timelike congruence that the density and pressure 
automatically satisfy a barotropic equation of state. In section 5, anisotropic fluid models are considered. 
In section 6, the two special physical subcases of a viscous fluid with heat conduction (satisfying a set of 
phenomenological equations) and a perfect fluid plus a magnetic field are investigated. In section 7, two 
component fluid models are considered, and particular attention is paid to the special physical subcases 
of a perfect fluid with pure radiation and two non-collinear perfect fluids (satisfying a variety of physical 
conditions such as, for example, separate energy conservation and linear equations of state). Finally, a single 
'tilting' perfect fluid is considered. We conclude with a discussion. 



2. General Equations 

Einstein's field equations are f] 



Gab = T ao (2.1) 



where the stress-energy tensor can formally be decomposed with respect to a timelike vector field u a ac- 
cording to 

Tab = ^u a u b + phab + q a u b + q b u a + ir ab , (2.2) 

where 

q a u a = 0, ir a a = 0, ir ab u b = , (2.3) 

and the projection tensor is defined by h ab — g a b + u a u b . In this formal decomposition /i, p, q a and ir ab are 
given by 

/i = T ab u a u b , (2.4) 
P = \h ab T ab , (2.5) 
q a = -h a c T cd u d , (2.6) 
TTab = h a c h b d T cd -^(h cd T cd )h ab . (2.7) 

For a fluid with four-velocity u a , these quantities denote the energy density, pressure, heat conduction and 
anisotropic stress, respectively, as measured by an observer moving with the fluid. 

In this paper, we consider spacetimes in which the shear, vorticity and acceleration (see jlj for 
definitions) of u a are all zero, whence the covariant derivative of u a can be written as 

u a -b = -9 Kb , (2.8) 

where 9 = u a . a is the expansion. The relevant equations in the case of zero shear, zero vorticity and zero 
acceleration are |l[] : 

the conservation equations 

A + (ji + p)0 + q a -,a = , (2.9) 

*We shall follow the notation and conventions in Ellis in particular, all kinematical quantities are defined therein. Also, 
Latin indices range from to 3 and Greek indices from 1 to 3, and subscripts indicate differentiation with respect to the 
relevant spacetime coordinates. 
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the Raychaudhuri equation 

the propagation equation 
the constraint equations 



h a b (p,b + TTb C ;c + % + 3 ft) = J (2.10) 

+ \o 2 + \ (p + 3p) = ; (2.11) 

n ab = 2E ab ; (2.12) 



^0, b h b a = q a , (2.13) 
H ab = ; (2.14) 

and the Bianchi identities (using the above) 

Ktb' b = \ht b V. b - ~0q t , (2.15) 

2 1 1 

+ 3^i"tro = -^h t a h m c q {a . c) + -q a - a h mt . (2.16) 

The tensors i? a 6 and i? a & are the electric and the magnetic parts of the Weyl tensor, C a bc d, and are given 

by 

E ab = C acbd u c u d , (2.17) 
H ab = \c acst u c r 1 st bd u d . (2.18) 

In the above, an overdot denotes differentiation along the fluid flow lines; for example, 6 = # a u a . 

The conditions of zero acceleration and zero vorticity imply that it is possible to choose a comoving 
coordinate system such that u a = 5® and such that the metric may be written as 

ds 2 = -dt 2 + g af3 {t,x"<)dx a dx . (2.19) 

The shear-free condition then implies that 

g af} = H 2 (t,x^)h a/3 (x^) , (2.20) 

and hence the expansion is 8 = 3H/H. We note that coordinates can also be chosen so that h a p is diagonal. 
Also, the Gauss-Codazzi equations imply that the Ricci tensor of the three-space of constant t is 

K/3 = *c0 + \ g a p{2^ - |fl 2 ) (2.21) 

which yields 

R* = 2(i - ^ 6 2 . (2.22) 

In the coordinates (t,x a ) of ( |2~19| ), u a = (1,0,0,0) , q = , 7r 0a = , 6 = d t 0, and h a b i>. b = 
implies ijj, a = 0. We shall adopt the coordinates (t, x a ) for the remainder of the paper. The energy density, 
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the pressure, the heat flux and the anisotropic stress tensor are then given by 



(2.23) 



P 
q a 



6 H 



H 



2(H/H), a = 2d t [d a (lnH)} = 2d a [d t (hxH)} 
1 



R. 



a(3 ~ g 9aj3 R* 



(2.24) 

(2.25) 
(2.26) 



Equations ( 2.23| ) and (2.26) follows directly from (2.21) and (2.22). Equation ( p. 24 ) follows by inserting 
(2.23) into the Raychaudhuri equation (2.11). The conservation equations ( fOi| ) and ( |2. 10 ) become 



A + -jf in + p) + g al3 q a ;p 

3H 



, 



dt(q a ) + ^ Vet + P,a + ff /37 7I"a/3; 7 = . 
Id 



(2.27) 
(2.28) 



Now, the three-metric g a p is conformal to the metric h a p , g a p — H 2 h a /3 , and thus the Ricci tensor R* a a is 
related to the Ricci tensor 3 R a /3, of the metric h a p Q by 



(2.29) 



= 3 R a p - V a V (laH) + Va(lnH)Vp(hiH) 

-[V 2 (lniJ) + V(hxff) • WQxiH)]h a p , (2.30) 



where V 2 iJ = h af3 V Q V H , Vif • Vi? = h al3 \7 a H\7 H h a ^ha 1 = 8%, and V Q is the covariant derivative 
with respect to the metric h a p- Thus, re-expressing (2.23)— (2.26 ) in terms of the three-metric h a p, we find 



P 

""a/3 



and 



where V" = h aP V f 



-H- 2 ( 3 R) - 2H- 3 V 2 H + H- 4 VH-VH + 'iH~ 2 H 2 

-2B- X B - X -ix , 
2d t [V a (lnH)] , 

3 R af 3 - \h af3 3 R - VaV/jQnff) + V a (\nH)W p {\nH) 
+ i [V 2 (lntf) - V(lnH) • V(lnff)] h a p , 



(i + SH^Hin + p) + H- l W a q a = , 
d t q a + 3H- 1 Hq a + V aP + H- 2 V f) TT af ] = , 



(2.31) 

(2.32) 
(2.33) 



(2.34) 



(2.35) 
(2.36) 



Finally, the Bianchi identities ( 2.15 ) and ( 2.16 ) reduce to 



H^V^n^ + H- s {W p H)-K a0 = g VaA* - H^Hq a , 
2d t n al 3 = -X7 a q p + 2q {a Vp ) {]nH) + - H 2 V 7 (if ~ 2 q 7 ) h a p 



(2.37) 
(2.38) 
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Equation (2.37) is equivalent to the contracted Bianchi identity V a {^R a p) — iV ) g( 3 i?), and ( 2.38Q simply 
expresses the result that dt(^R a p — \h a p 3 R) = 0. 



This completes our general analysis and any further progress can only be made on a case by case 



basis. 



3. Zero Anisotropic Stress 



I n th is secti on, we take the anisotropic stress tensor to be zero, n a b = 0. The equations (2.9)-( 2T6| ) and 
( [2.2l| )-( |2~22l ) then become 



A + O + PW + 9% = , 

4 „ 



p, a + q a + g v q 



o 



\o 2 + \ (At + 3p) = , 







E a b 

2 „ 
3 

2 (9«;/3 + <?/3;a) = ~ 9 7 ; 7 5a/3 > 



H a b — 
a Qa i 

= 0q a , 



i?* = 2/i — — I 



(3.1) 
(3.2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 

(3.8) 
(3.9) 



Equation (3.4) implies that the Weyl tensor is zero and thus the spacetime is conformally flat. Equation 
(3.6) can be integrated with the aid of (3.5) to obtain 



where / is an arbitrary function of integration. Comparing ( |3.1C| ) with (3.9), we find 

R* = f(t) . 



(3.10) 



(3.11) 



Therefore, each three-space (t = const) is a space of constant curvature and thus, we may choose coordinates 
so that the metric has the form ||, ||, 0] 



ds 2 = -dt 2 + fT 2 (dx 2 + dy 2 + dz 2 ) , 
VL{t,x a ) = a(t) {x 2 + y 2 + z 2 ) + b a (t) x a + c(t) , 

where a, b a and c are arbitrary functions of t. The Ricci scalar is then 

R* = 6(4ac- 8 afi b a bp) = R*(t) , 



(3-12) 
(3.13) 



(3-14) 
(3.15) 



[V is the covariant derivative with respect to the three-metric h a j3(— # a $)]. The expansion is now given as 



(3.16) 
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Equation (3.7) reduces to f2 >a/ a = |(V 2 f2)5 a/ 3 which is trivially satisfied by O as given in ( flf ). Also Q 
follows directly from (3.9) and the spatial gradient of (3.3). 



The density, /i, the pressure, p, and the heat flux, q ai can all be expressed in terms of f2 and its 
derivatives: 



pt = 2ft V 2 - 3Vft ■ vn 



ft 





2 U - 3 U -^v 2 o + v^.v^ 




The conservation equation ( |3.l| ) becomes 



A - 3(M + P)^ - ^ 2 (^ 2 i?* 



(3.17) 
(3.18) 
(3.19) 

(3.20) 



which is automatically satisfied by /i and p as given in (3.17) and (3.1 



The spherically symmetric solutions are found by setting b a (t) equal to zero in (3.13) and were 
originally discussed by Maiti 0. The metric can then be written as 



ds 2 = -dt 2 



dr 2 



[F{t)r 2 + G(t)] 2 



(3.21) 



where F and G are arbitrary functions. The model can then be interpreted as a perfect fluid with heat 
conduction, where the fluid four- velocity is u a = 81 . The density, pressure, and heat flux are now given by 



3(Fr 



G) 2 (Fr 2 + G)- 2 + Y1FG 



p = 2(Fr 2 + G)(Fr 2 + G)" 1 - 5(Fr 2 + G) 2 (Fr 2 + G)~ 2 - 4FG 



-4r[FG - GF](Fr 2 + G)' 



92 = <73 = 



(3.22) 
(3.23) 
(3.24) 



where qi is the heat conduction in the radial direction. We note that the model is conformally flat. Cos- 
mological models of this type have been studied by Kolassis et al. J9| and Banerjee et al. 

Perfect Fluid 

If bot h the anisotropic stress tensor and the heat flux vector are zero (7r a f, = q a = 0) then the stress- 
energy tensor (2^2) is that of a perfect fluid and the metric is the FRW metric (!]]. Thus if the stress-energy 
tensor is a perfect fluid, whose flow lines form a SIG timelike congruence, then the solutions of the Einstein 
field equations must be the FRW models. Furthermore, equations (3.2) and (3.6) imply that both p, and p 
are functions of the single variable i, [that is, p = ^{t)p = p{t)\ . Hence fi and p satisfy a barotropic equation 
of state p — p(p-). Therefore, p = p(fi) is a consequence of the assumptions and need not be specified a 
priori. 



4. Zero Heat Flux 



Here, we shall take the heat flux to be zero, q a = 0. We immediately obtain from equation (2.13) that 



9{t), and hence we can set H = H(t) without loss of generality. The anisotropic stress tensor is now 
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directly related to the Ricci tensor of the three- metric h a p\ equation ( [2.34 ) implies that 



TTa/3 



Raf3 — - h a /3 R 



The energy density and the isotropic pressure, as given by ( 2.3l| ) and ( 2.32| ), are 



H = 3H- 2 H 2 + -H- 2 ( 3 R) 



-IB- 1 H - H- 2 H 2 - 



'R 



The conservation equations 



and the Bianchi identities 



fi + 3H- 1 H{n+p) = , 

V Q p + H- 2 ^7T a0 = , 



1 



H- 2 V p n a p 



(4.1) 

(4.2) 
(4.3) 



(4.4) 
(4.5) 



(4.6) 
(4.7) 



are automatically satisfied. 



The zero heat flux models, under consideration in this section, can be subdivided into three distinct 
classes depending on the number of distinct eigenvalues of ir a p. If n a p has three distinct eigenvalues then 
the spacetime is of Petrov type I |n| . If TT al s has exactly two distinct eigenvalues then the spacetime is of 
Petrov type D ||o| [if] . Finally, if 7r Q ^ has only one eigenvalue (the eigenvalue is zero since ir aj 3 is tracefree, 
and hence ir a p is identically zero) then the spacetime is conformally flat and corresponds to a perfect fluid 
FRW model. 



In comoving coordinates, the metric has the following form: 

ds 2 = -dt 2 + H 2 (t)h af3 (x~ l )dx a dx p 



(4.8) 



Thus any spatial three-metric will give rise to a shear-free, irrotational, geodesic model with zero heat flux. 
In particular, there exist examples of spatially homogeneous spacetimes of all the Bianchi types [Q, I-IX, 
and the Kantowski-Sachs solution Jl3[ ] that admit a SIG timelike congruence. However, only Bianchi types 
I, V, VIIo, VII/j and IX can give rise to FRW models. Indeed, Mimoso and Crawford [jUj recently drew 
attention to the same fact: "No anisotropic model can simultaneously exhibit a perfect fluid matter content 
and a shear-free timelike congruence . . . But the possibility still exists of considering an imperfect fluid for 
shear-free anisotropic models" (where presumably the authors mean a shear-free timelikc fluid congruence, 
otherwise the statement is untrue — see section 6). 

As an illustrative example, we shall consider the following Bianchi IV spacetime with "equal scale 
factors," since this does not contain any FRW models as special cases: 



ds 2 = -dt 2 + H 2 (t){dx 2 + e 2x [x 2 + l)dy 2 + 2xe 2x dydz + e 2x dz 2 } 



(4.9) 



Taking the source to be an anisotropic fluid with the four- velocity of this fluid given by u a = 8ft, we 
observe that the acceleration, shear, and vorticity of the fluid congruence all vanish, and that the heat flux 
as measured by a comoving observer is zero. Using (4.1), we find that the non-zero components of the 
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anisotropic stress tensor are 



(4.10) 



' yy 



1 

3 ' 

J*(l 3 ?-2x-\\ , (4.11) 



3 3 



e 2x [^x~l) , (4.12) 

n zz = \e 2x . (4.13) 

The eigenvalue of ir a p can be found by solving the characteristic equation det(\h a p — ^aft) = 0. There are 
three distinct eigenvalues; namely, 

A = |, X=\±\Vl. (4.14) 



Thus, the spacetime is of Petrov type I. The energy density and the isotropic pressure are given by (4.2) 



and (4.3), respectively, with R = — 13/2. 



We note that for each three- metric h a p there exists an anisotropic stress tensor 7r Q( a, but not all of 
the Tr a f3 are related to physical relevant matter. Therefore, it is of interest to study models with some extra 
physical conditions. Perhaps the most reasonable condition to assume here is an equation of state between 
the energy density, /i, and the isotropic pressure, p. 

Equation of state: p = p(fi) 

Here, we consider the case where the energy density and the isotropic pressure satisfy an equation of state, 
p = p(fJ>), where both p and /i are differentiable. This implies that the integrability condition 

(V Q p)/i = (VcAOp ( 4 - 15 ) 
must be satisfied. Thus, inserting Q4.2| ) and ( |4.3| ) into ( 4.15| ) yields 



(V^dtOff- 1 ^) = . (4.16) 
Two possible solutions exist: either (i) \7 a [i = or (ii) dtiH^ 1 !}) = 0. 
(i) V Q /i = 

Here [i = n(t), and ((O) implies that the Ricci scalar is constant, 3 R =const. The isotropic pressure is 



now also a function of t, p = p(t). The anisotropic stress-tensor, as given by (4.1), is now subject to the 
constraint V a n a p = 0. Therefore, models with zero heat flux where the energy density and the isotropic 
pressure satisfy a barotropic equation of state are characterised by a single scale factor, H(t), and a spatial 
three-metric whose Ricci scalar is constant (not to be confused with three-spaces of constant curvature). 

We now make a brief comment on the consequences of the various energy conditions for the above 
models. The weak energy condition [|~5| — T a i,v a v b is non-negative for all unit timelike vectors v a — implies 
that the energy density, /i, is non-negative and hence, if 3 R is negative then (4.2) implies that either H > 



or H < for all time, t; thus, in the former case, the cosmology will be an open model. The strong energy 
condition — 2T a i,v a v b + T a a is non-negative for all unit timelike vectors v a — implies that [i + 3p > 0, 
and thus the deceleration parameters, q = —HH/H 2 , is non-negative. 
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(ii) d t (H~ 1 H)=Q 



Therefore, H = XH where A is a constant. We can integrate this equation to get 

H 2 = \H 2 + v , (4.17) 
where v is a constant. The energy density and the isotropic pressure are then 

H = - H- 2 ( 3 R + 6u) + 3A , (4.18) 
p = ~n - 2A . (4.19) 
If the fluid satisfies the strong energy condition then we can set A = —a 2 < and v — a 2 {3 2 > for some 



constants a and j3. If a ^ then (4.17) integrates to yield 

H(t) = /3sin(at) (4.20) 

which is a closed model. If A = then 

H(t) = jt (7 = const) (4.21) 
which is an open model. (In both cases, = —\-) 

Therefore, if the heat flux of the fluid is zero and the energy density and the isotropic pressure satisfy 
a barotropic equation of state then the hypersurfaces t = constant are for the most part hypersurfaces of 
constant Ricci scalar with the single exception of the equation of state ^2 = — i. 

Orthogonal Spatially Homogeneous Models 



The function H appearing in the metric ( [2.20 ) must be independent of the spatial coordinates for 



the metric to represent a spatially homogeneous model. Therefore, the expansion also only depends on t, 
9 = 0(t), and hence the heat flux is zero. 

To make further progress, we assume that the matter associated with these spatially homogeneous 
models inherits the symmetries [||; in particular, all scalars should be independent of the spatial coordinates 



Hence, both the energy density, fi, and the isotropic pressure, p, are now functions of only t. Equation (4.2) 
then implies that the Ricci scalar, 3 R, is constant, which is case (i) for zero heat flux models with a barotropic 
equation of state. The eigenvalues of the anisotropic stress tensor, ir a p, should also be constants. Indeed, 
as might be expected, all of the Bianchi types (I— IX) and the Kantowski-Sach models admit solutions with 



constant Ricci scalar such that the eigenvalues of Tr a p, as given by (4.1), are constant 



5. Anisotropic Fluid 



The stress-energy tensor for an anisotropic fluid is 

Tab = fJ,u a Ub + p\\ n a n b + p± (u a u b - n a n b + g ab ) , (5.1) 

where u a is a unit timelike vector, u a u a = —I, and n a is a unit spacelike vector, n a n a = 1, orthogonal to u a , 
u a n a — 0. The scalars pu and p± are the pressure parallel and perpendicular to n a , respectively. Equation 
(2.6) implies that the heat flux in the case of an anisotropic fluid is zero. Thus, shear-free, irrotational, 
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geodesic anisotropic fluids are a subcase of the models discussed in section 4. Equations ( p^ ) and (2.7) 
yield expressions for the isotropic pressure and the anisotropic stress tensor; 

P = \ Oil + 2 P^) , (5-2) 

TTab = (p\\ - P±){n a n b - ^(g ab + u a u b )} . (5.3) 



In terms of the comoving coordinates of ( 2.19 ), the spacelike vector n a — (0, n a ) and the anisotropic stress 
tensor is given by 

7Tq/3 = ip\\ -p±)(n a np - -H 2 h af3 ) . (5.4) 
The anisotropic stress tensor (^J) may be decomposed as |l(| 

7Ta/3 = P (N a N - ~ hap) , (5.5) 

where N a (x^) [= i2" -1 (*) 

n a ] is a unit vector with respect to the metric h a p and P is independent of t, 
P = P(x a ). The function P is a measure of the anisotropy of the fluid, 

P = H 2 (t) 0„ - p ± ) , (5.6) 

and is subject to the constraint equation 

iv Q ( 3 i?) = V p (PN p )N a - iv Q P + PN^VpNa . (5.7) 



The energy density and the isotropic pressure are given by expressions (4.2) and (4.3), respectively. The 
anisotropic pressures are related to p and P by 

P 

P± = P ~ ~MP ' (5 ' 8) 
2P 

P\\ = P+ JJp ■ ( 5 - 9 ) 
The anisotropic stress tensor (|5.5|) has only two distinct eigenvalues, |P and —iP (double), and 



hence the spacetime is of Pctrov type D (if P ^= 0). However, if P is identically zero then (5.6) implies 



that pu = p± and thus, the spacetime is then a perfect fluid FRW model. Conversely, if the metric is of 



the form ( |4.S| ) and 7r Q( g has exactly two distinct eigenvalues then Tr a p must necessarily be of the form (5.5) 
and hence the stress-energy tensor will formally be that of an anisotropic fluid with four-velocity u a = 6? 
p6| . Furthermore, if p and /i satisfy a barotropic equation of state, p = p(fi), then either the Ricci scalar, 
^R, is constant or |^ = _I ( see section 4). The case where the Ricci scalar, 3 R, is constant is the most 
interesting; in particular, there exists examples of spatially homogeneous anisotropic fluid models of Bianchi 
types II, III, VI and VIII, and a Kantowski-Sachs solution [pi. 



6. Special Cases 

Viscous Fluid 

For a viscous fluid with heat conduction the following phenomenological equations are satisfied: 

7T Qb = -2i]<7 ab = , (6.1) 

q a = ~ X h a b (T. b + Tu b ) = - X h a b T. b , (6.2) 

V = Pt - CO , (6.3) 
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where £ is the bulk viscosity coefficient (77 is the shear viscosity coefficient), p t is the thermodynamic pressure, 
X is the heat conductivity, and T is the temperature. These quantities are restricted by equations of state of 
the form C = C(p,T) 

[rj = Tj{fJL,T)], x = xb 1 ^) , p t = pt(/J,,T), and the various energy conditions. In addition, the Gibb's 
relation, the baryon conservation law and the second law of thermodynamics must be satisfied. 

Since 7r a & = 0, these models are a subcase of the models discussed in section 3, and hence there 



exists c oordi nates in which the metric is given by equations ( 3.12 ) and ( 3.13 ), and ju,, p and q a are given by 

(ra-ra. 



As an illustration, let us consider the case in which T = T(fi), whence 
Q = C(p-) ; X — x(m) an d Pt = Pt{n)- For example, at the high temperatures in the early universe when the 
energy density was dominated by relativistic species (that is, when the universe was radiation dominated), 
/i = ^gaT 4 = fi(T), where g = g(T) is th e nu mber of e ffective degrees of freedom contributing to the 
universe at temperature T. From equations (3.2) and (3.6), we obtain 



1 



whence 



-XT C 



dT 
^ dfi 



dT 
dfi 



= 0(p) 



(6.4) 



(6.5) 



since fi. a (q a 7^ 0), else the sp aceti me is a perf ect fluid FRW spacetime. This implies that f(t) 

constant in (3.10) and p = p(jj) from ( |6.3|) . Equation (3.3) then implies that fi is also a function of /i. 



Now, the right hand-side of the Gibb's relation, ^d (^) + (^) (where n is the baryon number 
density), is a perfect differential, which implies that either n = n(/i) (whence, from the baryon conservation 
law, n = —n9 , n and hence /i is a function of \i only), or we have the integrability condition 



whence from (|6.5|) we obtain 



rp dpt , 
T 1H = ^ 



P> + Pt 



XT 



. dT 
Pt) d^ 



dpt 
dfi 



(6.6) 



(6.7) 



Assuming that ^ > (and noting that (fx + Pt){xT) 1 > if the weak energy condition is satisfied) we 



thus have < 0. Otherwise (6.7) can be regarded as an expression defining x- 



Summarising, we have that 

6 = /i(m) , 
U = h(p) , 
[in addition to p = fs(fJ.)]- These imply the conditions 

dfi _ dh_dji_ 

cV* dt ~ dt 8x a ' *- l >w 



(6.8) 
(6.9) 



(6.10) 



which imply conditions on the (arbitrary functions in the) metric using ( |3.16| ), ( |3.17D and 

(pa. 



For example, in the case of spherical symmetry, when the functions b a (t) = in ( 3.13| ) whence the 
metric is given by (3.21), and [i,p and q, are given by equations (3.22), ( 3.23| ) and (3.24), respectively, ( |6.8| ) 
implies that 

F(t)G{t) = a , (6.11) 
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where a is a constant, whence 



_ 3(F 2 r 2 - a) 2 F 2 



Equations (|6.9[) and (6.10) (/i, r fj, t t — /it/x r ) then imply that 



and hence 



FF = F 2 



F = ce" f 



(6.12) 



(6.13) 



(6.14) 



where (3 and c integration constants. A constant rescaling of r can then be used to set c = a, whence the 
metric can be written as 

dr 2 + r 2 d9 2 + r 2 sm 2 9d(/) 2 

[ar 2 e^ + e"^] 2 ' (6 ' 15) 



ds z = -dt z 



Spacetimes of this form have been studied before |l7j. From equations ( ft.22| )-(3.24) we can then write 

1 „o 



/! = 12 a 



p = 2f3 2 - 4a - 



where 
and 



= 3/3(1 - ae 2 ^r 2 )(l + ae^'r 2 )" 1 
8a/3r 



Qi = 



[ar 2 e 131 + e' 131 ]' 



(6.16) 
(6.17) 

(6.18) 
(6.19) 



Perfect Fluid plus a Magnetic Field 

If the source of the gravitational field is a combination of a perfect fluid and a pure magnetic field then the 
stress-energy tensor has the form 



T ab = fiu a u b + p(g a b + u a u b ) + T ab 



MAG 



(6.20) 



If h a is the local magnetic field measured by u a (and the local electric field e a is zero) then the stress-energy 
tensor for a pure magnetic field is given by the Minkowski tensor |fl8f 



J-ab — A 



1 



(Jab - Ua U b ) h 



h n h h 



(6.21) 



where A is the magnetic permeability (assumed constant) and h 2 = h a h a (h a is orthogonal to u a , h a u a = 0). 
For instance, the above stress-energy tensor, ( |6.20| ), would be appropriate for a plasma in a strong magnetic 
field when the particle collision density is low ]19|, |20[ | . 

Equation (6.2C) is formally of the form of an anisotropic fluid pl| with energy density and isotropic 

pressure 



/' 

P 



P+\h 2 

6 



(6.22) 
(6.23) 
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and anisotropic stress tensor 



TT ab = \h z \ - {g a b + u a u b ) - n a n b 



where n a — h a /h. The anisotropic pressures are formally 



A 



P± 

P\\ 



P 



(6.24) 

(6.25) 
(6.26) 



The magnetic field must also satisfy Maxwell's equations which reduce to 

(u l h j - u j h% = 



(6.27) 



for a pure magnetic field |1S|]. If the integral curves of u 1 are shear-free, irrotational and geodesic then ( 6.27 ) 
implies 



h\_, 



, 



Contracting (6.29) with h % implies that 



- 9 hi + hi-j u 3 



' ' 2 Qh 2 = 



Now, equation (5.6) implies that 



2 W ' 3 



Yh 1 = P{x a )/H 2 {t) 



(6.28) 
(6.29) 

(6.30) 
(6.31) 



Inserting the result (3.31) into (6.30) yields H = 0, and thus the expansion is zero. Therefore, the magnetic 
field is static. In comoving coordinates, (4.8), equation ( 6.29| ) is trivially satisfied and ( 6.28| ) reduces to 



V a h a = . 

In addition, h a must also satisfy V^7r Q/ 3 = iVc/i, which reduces to 

7^- V a 7i = \h\7 a h — h^V ph a 



3A 



(6.32) 



(6.33) 



We note that there is no non-trivial solution if p is given by the equation of state p = (7 — with 7 
constant and in which the weak energy condition holds. 



7. Two Component Fluids 

Here, we consider the stress-energy tensor associated with a mixture of two fluids. We take one of the 
fluid components to be a perfect fluid whose flow forms a SIG timelike congruence (with fluid four-velocity 
u a = <5q); the density and pressure of this perfect fluid are denoted by fix and p\, respectively. We denote 
the contribution of this fluid to the stress-energy tensor by 



(1) T ab = (pi + pi)u a u b +pig a b 



(7.1) 
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The second fluid will be taken to be cither another perfect fluid with a velocity v a that is tilted with respect 
to u a or a pure radiation field with the following, respective, contributions to the stress-energy tensor 

(2) T ab = (fl 2 + P2)v a V b + P2 9ab , (7.2) 

™T ab = ek a k b . (7.3) 

where v a is a unit timelike vector (v a ^ u a ), and k a is a null vector, k a k a = 0. The total stress-energy 
tensor is simply 

T ab = (1) T ab + WToj . (7.4) 

The stress-energy tensor ( [7.4[ ) is equivalent to a single fluid with a non-zero heat flux and a non-zero 
anisotropic stress tensor. However, the anisotropic stress tensor associated with the fluid velocity u a is 
directly related to the heat flux (relative to u a ) via an equation of the form: 

7T Q b = TT^q a qb - ^{q c q c )(gab + w a u b )j . (7.5) 

In comoving coordinates (u a = 5§), the non-zero components of the anisotropic stress tensor are 

7i"a/3 = tt \q a qp - - {h lS q^q s ) h a[j 1 . (7.6) 



Letelier [22 has examined a two-perfect-fluid model of an anisotropic fluid. In particular, he 
examined the case where the stress-energy tensor consists of the sum of two perfect fluids, and of one 
perfect fluid and a null fluid. He studied the two-perfect-fluid model in the instance where both the perfect 
fluid components were irrotational. Ferrando, Morales and Portilla J23| have studied the two-perfect-fluid 
model where one of the fluid components was shear-free, irrotational and geodesic. They used an initial 
value formulation of general relativity to construct solutions. We essentially rederive their result using 
straightforward index notation, and indicate that the form of the metric they derive is appropriate for 
all shear-free, irrotational, geodesic timelike congruences when the anisotropic stress tensor has the form 
specified by equation (|7.5|). The results follow quite quickly from the general analysis of section 2. 



Ins erting the expression (7J3) for the anisotropic stress tensor into the left hand-side of (2.38) and 



then using ( p.3q ) to get an expression for dtq a yields 

= A a qp + Apq a + Bq a qp + C h a p (7.7) 

(the precise expressions for A a , B and C are unnecessary for the remainder of the discussion) . The unit 
spacelike vector 

w a ee ga (7.8) 

then satisfies 

V/3W a = -^@{h aj3 - w a w p ) + a a wp , (7.9) 



where 



6 = V a w a , (7.10) 
a a = w p Vf3W a [=A a -{h 01 A p w 1 )w a \ . (7.11) 

Therefore, the vector w a is shear-free and twist-free, in terms of the three-dimensional geometry. Thus, the 
three spaces admit an umbilical foliation p4| , and hence, there exists a coordinate system such that the 
three-metric has the following form jll], |25| 

h aP dx a dx fi = a 2 (x a )dx 2 + b 2 (x a )(dy 2 + dz 2 ) (7.12) 
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and the vector 



Equations (7.13) and (2.33) imply that H = H(t,x). Equation (2.34) then implies that 

3 R xy = V x V y {lnH) = -d x (lnH)d y (lna) , 
3 R XZ = V x V z (\nH) = -d x (\nH)d z (\na) . 



(7.13) 

(7.14) 
(7.15) 



But d x (hiH) is a function of both t and x (otherwise q a is identically zero). Therefore, we must have 
3 R xy = 3 R XZ = and hence, we can take a = 1 without loss of generality. Calculating 3 R xy and 3 R XZ for 



the metric (7.12) with a = 1, we find 

3 



R X y 
3 R X z 



(7.16) 
(7.17) 

(7.18) 

for some, as yet unknown, functions / and </>. Thus, the metric for the two component fluid (7.4) may be 
written as 



and hence 



= -b- 2 d xy b + b~ 3 d x bd v b 
= -b- 2 d xz b + b- 3 d x bd z b 

b = f(x)cj>(y,z) 



ds 1 



-dt 2 + H 2 (t,x){dx 2 + f 2 (x)(j> 2 (y,z) (dy 2 + dz 2 )} 



The only non-zero components of the anisotropic stress tensor, (2.34), for the above metric are 

Tt xx = -2f~ 2 (f)~ 2 iiyy = -2/~ 2 (j)~ 2 tt zz 

2 9 
= o ~ ~ H x — H~ H x f fx — f fxx + / " f x 



(7.19) 



(7.20) 



(7.21) 



We note that H is not a separable function, H ^ X(x)T(t), otherwise the energy flux and the anisotropic 
stress tensor would both be zero. 



Finally, ir xx must also satisfy 



7r q 



(7.22) 



where 

q = 2d t d x (lnH) (7.23) 

and the energy flux is simply q a — q5^. The functions q and it are therefore determined by (7.22) and ( |7.23| ). 
No further progress can be made in general, without specifying some sort of equation of state, usually an 
equation relating q and ir. For illustrative purposes, we will consider (i) a perfect fluid and a null fluid 
mixture, and (ii) a two-perfect-fluid mixture. 



Perfect fluid plus pure radiation 

The stress-energy tensor for a mixture of a perfect fluid and pure radiation is 

Tab = (fii + Pi)u a u b + pigab + ek a k b , 



(7.24) 



where k a is a null vector. We are free to rescale k a and e such that k a u a = — 1. Thus, in comoving 
coordinates we choose k a = (—1, k a ) and k a = H5^- Then, the various quantities (fi,p, q a , it a p) are 



V 
q a 

""a/3 



Mi 
Pi 
< //- 



1 

3 6 



eH 2 {S%Sl 



3 hap) 



(7.25) 
(7.26) 
(7.27) 
(7.28) 
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Therefore, 
and 



Using ( 7.29 ), we obtain 



q = eH = 2d t d x {\nH) 

2 2 
n xx = -eH 



k xx = -Hd t d x (\nH) 



Comparing ( |7.3l| ) to ( |7.2l| ) we find 



for some constant k, and thus 



+ 



i" 1 = i + l(y 2 + z 2 ) 



The functions H and / must also satisfy 

2Hd t (H- 1 H x ) + H- X H X 



2H- 2 H 2 X + H- 1 ^}- 1 ^ = 

I : f,r + f -7r - /••/ 1 



(7.29) 
(7.30) 

(7.31) 

(7.32) 
(7.33) 

(7.34) 



The above equation follows directly from (7.21) and ( [7.31 ). We seek solutions of (7.34) where H is not 
separable [H ^ X(x)T(t)]. It suffices to say that such solutions can be found for non-separable H. For 
example, suppose that / = 1 , k = and H = H(r) where r = x + at (a constant) then H(t) satisfies the 
differential equation 



d 2 H 2(<xff-l) fdH 



dr 2 H(2aH + 1) \dr 



Equation ( 7.35 ) can be integrated to yield 

H' 



dH\ (2aH + lf 



where c is an arbitrary constant. Thus, the line element 



ds 2 



dH 
IF 



■zH 2 



dx\ + H z (dx z + dy z + dz z ) , 



(7.35) 



(7.36) 



(7.37) 



where H' is given by ( 7.36| ) with a = 1, is an example of a spacetime where the stress-energy tensor consists 
of a perfect fluid (with 4- velocity u a — H'S^) plus pure radiation such that the perfect fluid's flow lines 
form a SIG timclikc congruence. 

Once (f>, f and H are known then the functions ni,pi, and e can be calculated from 

Hi = fi - e , (7.38) 
Pi = P ~ o e ' ( 7 - 39 ) 



= 2H- 1 d t (H- 1 H x ) , 



(7.40) 



where \i and p are given by (2.31) and (2.32). 
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Two-perfect-fluid mixture 



The stress-energy tensor for a two-perfect-fluid mixture is 

Tab = (pi + Pi)u a u b + (pi + p 2 )g a b + (A*2 + P2)v a v b 



(7.41) 



We shall assume that the two fluid four- velocities are not collinear, u a ^ v a , else the stress-energy tensor 
would formally be that of a single fluid with energy density /ii + ^2 and pressure p\ + p 2 and hence 
the resulting spacetime would be FRW. According to (2.4) and (2.5) the energy density and the isotropic 
pressure of the mixture are 



ju = Mi + (pa + p 2 ) cosh tp - p 2 
P = Pi 



where coshi/; 



P2 + o 1^2 + pz) sinh ip , 
-v a u a and ^ is called the tilt angle. The energy flux, ( [2.6| ), is 

9a = (P2 + P2) COs\ltp(v a + COSh^Ua) 



and the anisotropic stress tensor, (2.7), is 

1 



TTafc 



<7a qb 



^ (g cd qcqd)(gab + w Q "&) 



(7.42) 
(7.43) 

(7.44) 
(7.45) 



In comoving coordinates, u a 
(— cosh-0, sinh tp 5% ) , and 



cosh 2 ')/' (fJ.2 + P2) 

= —6®, whence the four- velocity of the tilting perfect fluid is given by v a = 
q a = (pi + Pa) cosh ip sinh ip 5* , (7.46) 



and 



1 



cosh %/}(ji2 + Pa) 



q a qp - \ [h lS q~f qs) h a p 



Therefore, the quantities q and ir are given by 

q = (A*2 



P2) cosh'0 sinhf/> 
1 



(/xa + P2) cosh ip 



From ( 7.48 ) and ( 7.49 ), we obtain the identity 



qn 



sinh ip 
cosh 



(7.47) 

(7.48) 
(7.49) 

(7.50) 



Now, the four quantities [i,p,q, and 7r are given, respectively, by equations (2.31), ( 2.32| ), ( 7.23| ) [see also 
(2.33)] and (7.22) [where tt xx is given by (7.21); see also (2.34)] in terms of the functions H,f and (j> [in 
metric (7.19)] and their derivatives. In addition, [i,p,q and 7r are given in terms of the five (unknown) 
physical quantities (J,i,pi, /xa,pa and tp through equations (7.42), ( 7.43j ), ( 7.48| ) and ( 7.49| ), respectively. At 
this point, therefore, the system of (four) equations is underdetermined, and no further progress can be 
made until physical conditions on fj,±,p±, /i2)Pa and tp are specified. If one such condition is specified, the 
remaining four physical quantities can then in principle be expressed on terms of H , f and <p and their 
derivatives. Two conditions on fii,pi, ^2,P2 an d tp would then [through equations (2.31), ( |2.32 ), (7.22) and 
( 7.23 )] give rise to a (differential) equation in terms of H, / and <p that would need to be satisfied (further 
restricting the form of the metric). 
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Let us consider the following conditions: (i) separate energy conservation, (ii) linear equations of 
state for p\ and P2, (hi) a constant tilt angle ip, (iv) the second perfect fluid is due to a scalar held, and (v) 
a single tilting perfect fluid. 

(i) Separate energy conservation. Here we assume that the energy momentum of each perfect fluid is 
separately conserved; that is, 



(A) T 



{A = 1, 2). Due to total stress-energy conservation, this leads to one extra constraint: 



Ml 



Mi + Pi 



(7.51) 



where p\ is a function of t only, p\ — pi(t). We note that if y,\ and p\ satisfy an equation of state of the 
form fix — fii(pi) — /ii(t), we immediately have that q a = 7r a & = and the resulting spacetime is FRW. 

(ii) Linear equations of state. Here we assume that each fluid satisfies a linear equation of state, viz., 

Pi = (7i , (7-52) 

P2 = (l2 - 1)M2 , (7.53) 

where the 71 and 72 are constants. Equations ( [7.42 ) and ( [7.43 ) then yield 

p + (1 - 71) M = M2{7i(72 - 1) + g72[-3 7i cosh 2 -!/' + 4 cosh 2 1/' - 1]} , 



and equations (7.49) and (7.50) imply 



1 



M2 = 



and 



1 



whence 



1 



72 



72 7r cosh ip cosh tp 

= 7T <^ p + (1 - 71) (Li - q 2 7T 



1 2 2 

= 1 — q 7T 



7i 



7i 
72 



(iii) Constant tilt angle. If if} = tpo = constant then from (7.50) we have that 

sinh -00 



q 7T 



cosh -00 

whence from ( 7.22| ) and (7.22) we have that 



o 



a const 



4a 



d t 8 x (lnH) , 



(7.54) 

(7.55) 
(7.56) 

(7.57) 
(7.58) 



which is similar to the expression given by (7.31) (with H 'replaced by a'). In particular, ix xx is a function of 
t and x only. Therefore, 4> again satisfies the differential equation (7.32) for a constant k, and is consequently 



given by equation (7.33) [that is, 



= 1 



j{y 2 + z 2 )], and now H and / satisfy an equation identical to 



(7.34) except that H is 'replaced by a' in the first term on the left hand-side [that is, the first term on the 
left hand-side of equation (7.34) is now 2ad t (H~ 1 H x )\. 

(iv) Perfect fluid plus a scalar field. The energy-momentum tensor of a scalar field is formally equivalent 
to that of a perfect fluid with 



M0 = \ (V^) 2 + V{<t>) , 



V(<f>) , 



(7.59) 
(7.60) 
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where V is the potential of the scalar field. Therefore, a perfect fluid plus scalar field source can be formally 
treated as a two-perfect-fluid mixture (with separate energy conservation). The field 4> satisfies a (Klein- 
Gordon) differential equation; therefore, for a specific V((f>) we will obtain a relation (effective 'equation of 
state') between and p$. 



(v) Single tilting perfect fluid. Here we assume that 

Hi = pi = 



(7.61) 



so that the energy-momentum tensor represents a single perfect fluid whose four-velocity is tilting with 
respect to the shear-free, irrotational and geodesic timelike congruence (and hence the results at the end of 
section 3 do not apply). 

From (|7.42[) and (|7.43|) , using ( |7.6l| ), we have that 



M + P = g (j"2 + P2) [3 cosh V + sinh tp] 



whence from ( |7.48| ) and fl7.49| ) we obtain 



and hence, we have the identity 



3(/it + p) cosh?/; sinh-0 
3 cosh 2 ip + sinh 2 ip 

3 cosh 2 tp + sinh 2 tp 
3(/k + p) cosh 2 tp 



irq 2 [(jjt + p) - — TT q 2 ] 



Now, from (2.31) and (2.32) we obtain [for metric ( |7.19| )] 



(7.62) 

(7.63) 
(7.64) 

(7.65) 



P = -2H- 1 H tt + 2H- 2 H 2 + -H- 2 {-2H- 1 H xx + H- 2 H 2 X 

-ah- 1 h x /- 1 f x - 2/- 1 f xx - r 2 f x + r 2 $( y , z )} , 



where 



"yy 



+ <!>„) + {4> 2 y + €) 



Using equations (7.2C)-(7.23) and (7.66), we see that equation ([7.65]) is of the form 



a(t, x) $ 2 (y, z) + b(t,x)®(y,z) + c(t,x) 







(7.66) 
(7.67) 

(7.68) 



and hence <&(?/, z) = k, a constant; consequently <j) again satisfies the differential equation (7.32) and is thus 



given by equation ( 7.33j ) [that is, 



= 1 + f (y 2 + z 2 )] 



Finally, using (fj\2|) and (^2|), fi\6|) becomes 



H x t 
If 



H x H t 
~H 2 ~ 



(P + P) 



(7.69) 



where (/i + p) is given by ( |7.66| ), tt xx is given by ( [7.20 ) and ( 7.21 ), and where $(y, z) — k. Equation ( 7.69| ) 
is a differential equation for H(t, x) and f{x). 

Similar to the case of a perfect fluid with pure radiation, it is possible to find special solutions to 
the differential equation ( |7.69| ). If / = 1, k = and H = H(t), where r = x + at (a constant), then ( 7.69| ) 
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becomes a quadratic equation in H" /{H') 2 (H' = ^p): 



H 2 {(1 + 6a 2 ) if 2 - 4} 



H" 

Wf 



+ 2H {(2 - lha 2 )H 2 - 3} 
+ 4 {(1 + 6a 2 )H 2 + 1} = 



H" 



(7.70) 



This quadratic equation can be solved and the solution for H(t) can be formally obtained. As an illustration, 

4 
3 ' 



for the specific choice of a 2 = 4 



H" 18H 2 + 3- V216i/ 2 + 25 



(ij') 



/\2 



H{9H 2 - 4) 



is one of the solutions of the quadratic equation (7.70). Equation (7.71) can be integrated to yield 



H' 



\9H 2 - 4|- |5 V216ff 2 + 25 - 25 - 108ff 2 | 



H 2 |11 V216i? 2 + 25 - 73 - 108 # 2 | 
where c is an arbitrary constant. Thus, the line element 

3 fdH 



ds 2 



4 V H' 



dx 



H 2 (dx 2 + dy 2 + dz 2 ) , 



(7.71) 



(7.72) 



(7.73) 



where H' is given by j7.72j ), is an example of a model with a single tilting perfect fluid where the fluid is 
tilting with respect to a shear-free, irrotational, and geodesic timelike congruence. 



Once , / and H are known, the scalar \ii and P2 can be calculated from 

3 



(l, 2 
P2 



A* 2 ^ xx ' 
2q 2 



37T., 



(7.74) 
(7.75) 



where fj,, tt xx and q are given by equations (2.31), (7.21) and (7.23|), respectively. 



8. Discussion 



In this paper, we have considered spacetimes, with a general energy-momentum tensor [formally decomposed 



with respect to u a according to (|2.2j)-(2.7)1, in which the shear, vorticity and acceleration of a timelike 
congruence u a are all zero. In such spacetimes, it is always possible to choose a synchronous coordinate 
system with u a = Sf [see ( 2.19| )— ( 2.20[) ] such that all dynamical information is encoded in a single scalar 



function H(t,x a ) and such that the geometry of each of the spacelike hypersurfaces t — const is encoded 
(up to a conformal factor) in the positive-definite three- metric ha^x 1 ). 

In section 3, spacetimes with zero anisotropic stress (with respect to u a ) were considered. Each 
three-space (of constant t) was shown to be a space of constant curvature. However, the curvature was not 
necessarily the same on all the t = const hypersurfaces. In section 4, spacetimes with zero heat flux vector 
(with respect to u a ) were considered. The scalar function H was demonstrated to be independent of the 
spatial coordinates [that is, H = H(t)], and the anisotropic stress tensor was shown to be given entirely 
in terms of h a p(x' y ) through ( |4.l| ). When the energy density and isotropic pressure satisfy a barotropic 
equation of state, p = p(/i), we saw that (except in the special case ^ = — |) the Ricci scalar of the 
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three-surface [with metric ha^x 1 )} is constant, whence the anisotropic stress-tensor is 'divergence-free' (all 
of the Bianchi and Kantowski-Sachs models admit subcases which fall into this category). In section 5, an 
anisotropic fluid source, (5.1), was studied. These models were shown to be a subset of the models discussed 
in section 4. The fluid's anisotropic stress tensor was found to possess two distinct eigenvalues (and hence 
the spacetime is of Petrov type D). 

In section 6, two physically relevant energy-momentum tensors were considered, which were sub- 
cases of models dealt with in earlier sections. First, a viscous fluid with heat conduction satisfying the 



phcnomcnological equations (6.1)— (6.3) was considered, a subcase of the models with zero anisotropic stress 
tensor. In the case that the temperature depends only on the energy density, T = T(ji), the various govern- 
ing equations and thermodynamical laws were found to imply a number of constraints; in particular, we were 



able to 'integrate' these conditions for the spherical symmetr ic cas e, ( 5.15)— (6.19). Second, a perfect fluid 
plus a magnetic field with energy-momentum tensor given by ( 6.20 ) and (|6 . 21 ) , which is formally equivalent 
to an anisotropic fluid, was considered. For this case, the Maxwell equations were shown to imply that the 
expansion must be zero and hence, the magnetic field must be static. 

Finally, in section 7, a number of physically relevant cases were investigated which can be considered 
(or grouped together) as two component fluids in which the stress-energy tensor is given by (7.4) with 
(7.1) (representing a perfect fluid whose flow forms a SIG timelike congruence) and one of (7.2) or (7.3) 
[representing a second (tilting) fluid] . The stress-energy tensor is formally equivalent to that of a single fluid 
with a non-zero heat flux and a non-zero anisotropic stress tensor that is directly related to the heat flux via 
(7.5). In this particular case, coordinates can be chosen such that the metric has the simplified form ( 7.19| ) 
and only depends on three arbitrary functions, H(t,x), f(x) and (f>(y,z). A number of special cases were 
then considered. First, the 'second fluid' was taken to be pure radiation, ( |7.3| ), so that the total stress-energy 
tensor is due to a mixture of a perfect fluid and pure radiation. In this case, the function <j)(y, z) is given 
by equation (7.33), and the functions H(t,x) and f(x) must satisfy the differential equation (7.34). In the 
second case, the source was taken to be two non-collinear perfect fluids, (7.41). The various quantities in 
the formal expression for the stress-energy tensor were then found to be given by equations ( [7.42 ), (7.42), 
( [7.48 ) and (7.49). These four equations relate five physical quantities to the metric functions H,f and 0, 
and their derivatives. Consequently, no further progress could be made until additional conditions were 
imposed on the five physical quantities. A number of conditions were then considered. The case of a single 
perfect fluid tilting with respect to the SIG timelike congruence is of special interest. In particular, we note 
that a general relativistic spacetime that admits a SIG timelike con gruen ce and whose stress-energy tensor 
is a (single) perfect fluid need not be FRW. The metric given by ( 7.73 ) represents a specific example of 
such a spacetime. However, the spacetime is an FRW model when the perfect fluid's flow lines form a SIG 
timelike congruence. 
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